Hole Localization in Underdoped Superconducting Cuprates Near l/8th Doping 
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Measurements of thermal conductivity versus temperature over a broad range of doping in 
YBa2Cu306+a; and HgBa2Ca n _iCu n C>2,i+2+,5 (ra=l,2,3) suggest that small domains of localized 
holes develop for hole concentrations near p=l/8. The data imply a mechanism for localization that 
is intrinsic to the CuCb-planes and is enhanced via pinning associated with oxygen- vacancy clusters. 

PACS numbers: 74.72.-h, 74.25.Fy, 74.25.Bt, 74.25.Dw 
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There has been considerable recent interest in novel 
charge- and spin-ordered phases which may compete with 
superconductivity in the cuprates JjJ . This stripe order is 
favored at CuC>2 planar hole concentrations near p=l/8 
for which the modulation wavelength is commensurate 
with the lattice, and when an appropriate pinning poten- 
tial is present. Inelastic neutron scattering results M sug- 
gest that stripe modulations are disordered and/or fluc- 
tuating in Lai_ x Sr x Cu04 (La-214) and YP^CuaOg+x 
(Y-123). Nuclear magnetic and quadrupole resonance 
(NMR and NQR) studies ||[| indicate the presence of 
localized holes in the CuC>2 planes. Whether localized 
holes are a general feature of cuprates and contribute 
to short-ranged charge/spin segregation and the normal- 
state pseudogap [^|^] are fundamental issues of current 
interest. 

Here we report the observation of anomalies at p—1/8 
revealed in the oxygen doping dependence of thermal con- 
ductivity (k) in Y-123 and Hg cuprates. Our principal 
finding is that a fraction of the doped holes in the CuC>2 
planes become localized near this particular doping level. 
This fraction correlates with the oxygen vacancy concen- 
tration in the charge reservoir layers, consistent with the 
formation of nanoscale, hole-localized domains that are 
pinned near oxygen-vacancy clusters. 

Two prominent characteristics of the in-plane heat con- 
ductivity in superconducting cuprates are a normal- 
state k predominantly of lattice origin, and an abrupt 
increase in k at T < T c . To motivate our proposal that 
k probes localized holes, we demonstrate for Y-123 that 
its magnitude and temperature derivative at T c corre- 
late with independent measures of local structural dis- 
tortions (planar Cu NQR) and the superfluid fraction 
(specific heat jump), respectively, throughout the under- 
doped regime (p < 0.16). Consider data from previous 
work H that illustrate the p=l/8 features in k (Fig. [I]). 
The upper panel shows normal-state data, tt(p, T=100K), 
for an Y-123 polycrystal and the ab-plane of single crys- 
tals. The lower panel shows, for the same specimens, 
the dimensionless slope change in k at T c defined as || 
T = -d{n s /K n )/dt\ t=1 , where t = T/T c and K s {n n ) is the 



thermal conductivity in the superconducting (normal) 
state, r measures the change in scattering of heat carri- 
ers (electrons and phonons) induced by superconductiv- 
ity. Both quantities exhibit local minima near p=l/8. 

Figure [l] (a) implies a sharply reduced lattice con- 
ductivity (kl) about p=l/8, since the electronic thermal 
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FIG. 1. (a) Doping dependence of the normal-state re at 
T=100K for Y-123 polycrystals and the ab-plane of single 
crystals from Ref. || . Values of p were determined from ther- 
mopower measurements following Ref. ji| Also plotted (x's, 
right ordinate) is (0.08 + 0.16/jvq.r) -1 , where Inqr is the 
anomalous planar 63 Cu NQR signal (Ref. jj, see text). The 
dashed curve is a guide to the eye. The solid curve is dis- 
cussed in the text, (b) The normalized slope change in re(T) 
at T c vs doping for each of the Y-123 specimens in (a). Also 
shown are the normalized electronic specific heat jump, A7/7, 
from Ref. [p| (x's), and the /iSR depolarization rate (in fis -1 ) 
from Ref. 112] (solid curve), divided by 1.4 and 2.7, respec- 
tively, and referred to the right ordinate. The dashed curve 
is described in the text. 
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tivity (re e ) at T > T c in optimally doped, single-crystalY- 
123 represents only ~1Q% of the total re = re e + k,l [ jio| . 
We attribute this suppression in to phonon scattering 
by local distortions of the CuO polyhedra. In support 
of this interpretation is the doping behavior of the rel- 
ative intensity {Inqr) of anomalous planar 63 Cu NQR 
peaks for Y-123, indicative of charge inhomogeneity in 
the planes, and attributed to the presence of localized 
holes The simple assumption that the total thermal 
resistivity is a sum of a term proportional to this disorder 
(oc Inqr), and a constant term from other scattering, re- 
produces k(p) quite well [crosses, Fig. |l| (a)]. Studies of re 
in lightly-doped insulating cuprates and structurally sim- 
ilar manganites fill provide further evidence that bond 
disorder associated with localized charge is an effective 
phonon scattering mechanism in perovskites. 

The electronic specific heat jump for Y-123 ||, 
A7/7Q?) [crosses, Fig. |l| (b)], correlates remarkably with 
r(p). Though the origin of the slope change in re at 
T c (i.e. electronic or phononic) has been a subject of 
some debate, this correlation demonstrates unambigu- 
ously that r, like A7/7, measures the difference between 
the normal- and superconducting-state low-energy elec- 
tronic spectral weight. This conclusion is robust since the 
correlation holds for both polycrystals and single crys- 
tals; in spite of additional thermal resistance due to c-axis 
heat flow and granularity in the polycrystal, changes with 
doping reflect the intrinsic behavior of the in-plane heat 
conduction. 

The muon spin rotation (/jSR) depolarization rate (do) 
]l2| |, a measure of the superfluid density, matches the 
data for T and A7/7 at p < 0.09 and p > 0.15 when suit- 
ably scaled [solid curve in Fig. g (b)]. The sharp rise m 
Co for p > 0.15 is due to the superconducting condensate 
on oxygen-filled chains |l2l . The range of p over which 
r and A7/7 deviate from the fj.SH curve coincides with 
the range of suppressed re. The suppressed transfer of 
spectral weight below T c is consistent with the localiza- 
tion of a fraction of planar holes. That oq is not also 
suppressed near p—1/8 indicates that hole- localized do- 
mains do not inhibit the formation of a flux lattice in ad- 
jacent regions where holes are itinerant. The scenario we 
have outlined resolves the prior inconsistency of the non- 
monotonic A7/7 doping behavior with the continuously 
rising pseudogap energy scale, inferred from numerous 
experiments in the underdoped regime ||p|. 

If our interpretation is correct the maximum density 
of localized holes occurs at p — 1/8. This suggests a pos- 
sible connection with the physics of stripe formation and 
raises questions about the generality of this phenomenon 
and the role of oxygen vacancies. To address these is- 
sues we have studied re in HgBa2Ca m _iCu m 02 m +2+5 
[Hg-1201 (to=1), Hg-1212 (m=2), Hg-1223 (m=3)]. A 
single HgOa layer per unit cell contributes charge to 777 
planes in Hg-12(777-1)777 so that the oxygen vacancy con- 
centration, 1 — <5, increases with decreasing m (at opti- 



mum doping, <5 G pt ~ 0.18 (Hg-1201), 0.35 (Hg-1212) and 
0.41 (Hg-1223) ||). 

The preparation of the polycrystalline starting ma- 
terials is described elsewhere [Q. Transport measure- 
ments were performed on specimens (~ 1 x 1 x 3mm 3 ) 
cut from as-prepared disks; each was measured repeat- 
edly up to ten times after successive heat treatments at 
300-350°C in flowing argon or vacuum to reduce the oxy- 
gen content. The thermopower (S) and thermal conduc- 
tivity were measured during each experimental run us- 
ing 25/im type-E thermocouples and a steady-state tech- 
nique. Heat losses via radiation and conduction through 
the leads (10-15% near 300K and 1-3% for T < 120K) 
were determined in separate experiments and the re data 
corrected. Uncertainty in the contact geometry intro- 
duces a ±5% inaccuracy in re. 

Figure | shows S{T) and re(T) for the Hg-1212 spec- 
imen. The thermopower curves are labeled by p, as 
determined from S'(290K) and the universal S(p) re- 
lations established by Tallon et al ||. Data for Hg- 
1201 and Hg-1223 yield similar sets of curves. Fig- 
ure || shows T c (p)/T™ ax for each of the materials 
studied. The data follow the generic high-T c phase 
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FIG. 2. (a) Thermoelectric power and (b) thermal con- 
ductivity vs temperature for Hg-1212. Data at p = 0.121 
and p = 0.132 are omitted from both (a) and (b) for clarity. 
The inset shows re(T) near T c and extrapolated re™ curves for 
p = 0.153. 



2 




FIG. 3. T c (p) normalized to T™ ax for each of the Hg com- 
pounds. The solid curve is T c /T™ ax = 1 - 82.6(p - 0.15) 2 . 

curve j^j, approximated well by the parabolic form, 
T c (p)/T c max = 1 - 82.6(p - 0.15) 2 (solid curve, Fig. §). 

The doping dependencies, n(p 7 T '/T"^ x = 1.2) and 
T(p) for each material are shown in Fig. [| (a) and (b), re- 
spectively. As for the Y-123 data [||, k" was determined 
by polynomial fits to the normal-state data at t > 1.2 
(see inset, Fig. ||). Uncertainties in the heat-loss correc- 
tions and extrapolations of K n in the range 0.9 < t < 1.2 
are reflected in the error liars for T. Uncertainties in the 
extrapolations were established by varying the range of 
data and polynomial order used for fitting. 

The data for the Hg materials are qualitatively similar 
to that of Y-123, but with the magnitude of suppres- 
sion near p=l/8 in k and T increasing with the oxygen 
vacancy concentration. Consider first the T(p) data for 
Hg-1223 which falls on an inverted parabola, centered 
near p — 0.16 [solid curve, Fig. [| (b)]. This behavior is 
similar to the /iSR curve for Y-123, and suggests that 
the superconducting condensate is not suppressed in Hg- 
1223. An extension of our interpretation for Y-123 im- 
plies that very few holes are localized in this compound. 
The T data for underdoped Hg-1201 and Hg-1212 are 
suppressed relative to that of Hg-1223 and maximally so 
near p—1/8, consistent with the localization of a fraction 
of holes. 

A similar behavior is seen in the k data. For p < 0.16, 
Hg-1212 has n(p) indistinguishable from that of Hg-1223, 
except in the narrow range of p near 1/8 where it is 
sharply suppressed, k for Hg-1201 is substantially sup- 
pressed over a broader doping range. The implication 
is that in the absence of localized holes, k for under- 
doped Hg-1201 and Hg-1212 would follow the solid curve 
in Fig. U (b), an interpolation through a composite of 
data for all three compounds, excluding the data point 
at p ~ 0.18 for Hg-1223 and ranges of p for Hg-1212 and 
Hg-1201 where n is suppressed. The increase in n(p) in 
the overdoped regime (p > 0.16) for all compounds is 
presumably due to the rising n e , and possibly to a ris- 
ing Ki as well, the latter due, e.g., to a decrease in the 



phonon-electron scattering. 

Our data suggest that the maximum in n(p) observed 
near p=l/8 for Hg-1223 is an intrinsic feature of kl 
in underdoped cuprates. From kinetic theory, kl = 
1/3Clv 2 t, where Cl is the lattice specific heat, v is the 
sound velocity, and r the phonon relaxation time. Inter- 
estingly the behavior of the normal-state elastic constants 
(oc v 2 ) of single-crystal La-214 |lj] parallels that of k for 
Hg-1223. A substantial hardening of the lattice in the un- 
derdoped regime, with a maximum near p — x —1/8, is 
observed for all symmetries (shear mode data are shown 
in the inset, Fig. ||). An increasing v 2 is expected quite 
generally from elastic coupling to conduction electrons 
for a decreasing electronic density of states (the pseu- 
dogap) [jLfSj. The appearance of a maximum suggests 
a competing effect, possibly related to the increasingly 
magnetic character of the planes. Thus the maximum 
in k for underdoped Hg-1223 is plausibly attributed to 
a renormalized dispersion, and its suppression in com- 
pounds with localized holes to increased scattering. In 
spite of the oversimplification implicit in kinetic theory, 
these observations suggest that in the absence of local- 
ized holes, k(p) for Y-123 would have a similar behavior 
[solid curve, Fig. [I] (a)]. 

Our interpretation dictates a proportionality be- 
tween the suppression of n and T, allowing for 
a self-consistency check. The lattice thermal re- 
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FIG. 4. (a) Doping dependence of the normal-state re at 
rp jrjvmax _ 2 for each of the Hg compounds. Dashed curves 
are guides to the eye and the solid curve is described in the 
text, (b) The normalized slope change in re(T) at T c vs doping 
for each of the Hg compounds at the same doping levels as in 
(a). The solid line is 1.71- 250(p-.157) 2 . Dashed curves are 
described in the text. The inset shows elastic constant data 
from Ref. g. 
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sistivity (W ~ 1/k) introduced by randomly-distributed, 
hole-localized domains should be proportional to their 
volume fraction since k is a bulk probe. We then ex- 
pect r (and A7/7) to be suppressed by an amount pro- 
portional to the same fraction, i.e. AT(p) — aAW(p), 
where Ar and AW are the amounts by which the mea- 
sured r and W differ from values given by the solid 
curves in Fig.'s [l] and |I| This simple relation agrees 
very well with the data [dashed curves in Fig.'s |l| (b) 
and I (b)] using a 1212 ~ 1.7W/mK, a 1201 ~ 2.6W/mK, 
and cP^ 123 ~ 2.8W/mK. 

The Hg data make it clear that oxygen vacancies 
play only a supporting role in the localization of holes. 
The importance of l/8th doping implies that the phe- 
nomenon involves an excitation of the CuC>2 planes that 
is commensurate with the lattice. A plausible candidate 
is a small domain of static stripe order Q, nucleated 
via pinning by a vacancy-induced mechanism. Recent 
Raman scattering studies of optimally-doped Hg com- 
pounds JL7| implicate vacancy clusters in pinning: the 
oscillator strength of the 590 cm -1 mode, attributed to 
c-axis vibrations of apical oxygen in an environment of 
four vacant nearest-neighbor dopant sites, scales with 
Ak and Ar. Oxygen vacancies in both Y-123 and Hg 
cuprates |ll|Jl8| displace neighboring Ba 2+ toward the 
CuC>2 planes. Resulting local distortions of the CuO 
polyhedra could pin a charge stripe as do the octahedral 
tilts in (La,Nd)-214 |Q. Alternatively, a magnetic mecha- 
nism is possible. In the Hg cuprates the shift of Ba atoms 
associated with a cluster of four vacancies nearest a CuO 
polyhedron will induce a positive potential in the planes 
that inhibits its occupation by a hole, thereby suppress- 
ing spin fluctuations and possibly fixing a Cu 2+ spin at 
the site. Larger vacancy clusters, surrounding adjacent 
CuO polyhedra oriented along (100) directions, may in- 
duce several spins to order antiferromagnetically via su- 
perexchange. The presence of this spin-chain fragment 
would favor charge/spin segregation that is characteris- 
tic of stripe order. 

Consider the relevant length scales. If randomly dis- 
tributed, static stripe domains are to scatter phonons, 
their in-plane extent must be less than the phonon mean- 
free-path, A, and their mean separation comparable to A. 
Expressing A -1 as a sum of terms for scattering by these 
domains and by all other processes, A -1 = A^ r + A~ t ^ er , 
we use the in-plane thermal resistivity for Y-123 and ki- 
netic theory to find @ A sir ~ [3/C L vAW(p=l/8)] ~ 
70A as an estimate of the separation between domains at 
p=l/8. We can make a rough estimate of the fractional 
area of the planes having static stripes by taking the typ- 
ical domain size to be 2a x 8a (a is the lattice constant), 
the stripe unit cell suggested M for (La,Nd)-214. This 
yields a fraction 16(a/A str ) 2 ~ 0.05. The T data suggest 



that this fraction is somewhat higher in Hg-1201, but 
apparently below the 2-D percolation threshold of 0.50. 
This presumably explains why no substantial T c suppres- 
sion is observed near p—1/8 in Y-123 or the Hg cuprates, 
in contrast to the case of (La,Nd)-124; stripe domains in 
the latter system have longer-range order (> 170A) §. 

The authors acknowledge helpful comments from J. 
Ashkenazi, S. E. Barnes, M. Peter, and F. Zuo. Work 
at the University of Miami was supported by NSF Grant 
No. DMR-9631236, and at the University of Houston by 
NSF Grant No. DMR-9500625 and the state of Texas 
through the Texas Center for Superconductivity. 



[1] 

[2] 

[3] 
[4] 

[5] 

[6] 



I'] 



J. M. Tranquada et al, Nature 375, 561 (1995); Phys. Rev. 
Lett. 78, 338 (1997). 

J. M. Tranquada, ANL preprints, cond-mat/970211'; and 



cond-mat/9709325 



523 



[8] 

[9] 
[10] 

[11] 



[12] 
[13] 



[14] 
[15] 

[16] 

[17] 
[18] 
[19] 



P. C. Hammel et al, Phys. Rev. B 57, R712 (1998). 
P. C. Hammel and D. J. Scalapino, Phil. Mag. 74, 

(1996) . 

J. W. Loram et al, Phys. Rev. Lett. 73, 1721 (1993); J. 
W. Loram et al, J. Supercond. 7, 243 (1994). 

B. Battlogg et al, Physica C 235-240, 130 (1994); B. N. 
Basov et al, Phys. Rev. B 50, 3511 (1994); J. L. Tallon 
et al, Phys. Rev. Lett. 75, 4414 (1995); A. Loeser et al, 
Science 273, 325 (1996); G. V. M. Williams et al, Phys. 
Rev. Lett. 78, 721 (1997). 

C. Uher, in Physical Properties of High Temperature Su- 
perconductors, edited by D. M. Ginsberg (World Scientific, 
Singapore, 1993), Vol. III., p. 159. 

C. P. Popoviciu and J. L. Cohn, Phys. Rev. B 55, 3155 

(1997) ; J. L. Cohn, ibid. 53, R2963 (1996). 

J. L. Tallon et ai, Phys. Rev. B 51, 12911 (1995). 
K. Krishana et ai, Phys. Rev. Lett. 75, 3529 (1995). 
J. L. Cohn, C.-K. Lowe-Ma, and T. A. Vanderah, Phys. 
Rev. B 52, R13134 (1995); J. L. Cohn et ai, Phys. Rev. B 
56, R8495 (1997). 

Y. J. Uemura et al, Phys. Rev. Lett. 66, 2665 (1991); J. 
L. Tallon et ai, ibid. 74, 1008 (1995). 

O. Chmaissem et al, Physica C 217, 265 (1993); E. V. An- 
tipov et al, ibid. 218, 348 (1993); Q. Huang et al, Phys. 
Rev. B 52, 462 (1995). 

Q. M. Lin et ai, Physica C 254, 207 (1996). 

S. Sakita et ai, Physica B 219-220, 216 (1996); M. Nohara 

et ai, Phys. Rev. B 52, 570 (1995). 

M. Peter and E. R. Walker, J. Appl. Phys. 48, 2820 (1977); 

B. Liithi, J. Mag. Mag. Mater. 52, 70 (1985). 

X. Zhou et al, Phys. Rev. B 54, 6137 (1996). 

J. D. Jorgensen et ai, Phys. Rev. B 41, 1863 (1990). 

We take Ci,(100K)=1.3x 10 6 J/m 3 K (Ref. g) and v = 3 x 

10 3 m/s estimated from dispersion curves [L. Pintschovius 

et ai, Physica C 185-189, 156 (1991)]. 



4 



